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We use the weakly nonlocal hydrodynamics approach to obtain a dynamical equation for the 
peculiar velocity field in which the viscosity term is physically motivated. Based on some proper- 
ties of the Ginzburg-Landau equation and the wave mechanics analog of hydrodynamics we find 
the nonlocal adhesion approximation taking into account the internal structures of the Zeldovich 
pancakes. If the internal structures correspond to significant mesoscopic fluctuations, viscosity is 
probably driven by a stochastic force and dynamics is given by the noisy Burgers equation. 
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II. HYDRODYNAMICS AND ADHESION 
APPROXIMATION 



I. INTRODUCTION 



Let us define the peculiar velocity field as the vector 
field u = r Hr (where H = a/a and a is the scale 
factor). Consider the fluid motion near a point 0. Then 
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The complexity of the nonlinear structure formation 
process is not possible to be studied just using analytical 
methods. But the comparison of analytical approxima- 
tions and N-body experiments may help us to under- 
stand the physics of highly nonlinear collective phenom- 
ena. The Zeldovich approximation was the first analyti- 
cal try to predict large-scale structures like pancakes, fil- 
aments and voids . The problem with this kinematical 
model is that particles pass through the pancakes instead 
of clustering into smaller objects like groups and galax- 
ies. An important improvement of the Zeldovich model 
is the adhesion approximation in which a viscous 
term is introduced to prevent particle crossing. However, 
this modification is not intended to describe the internal 
structure of pancakes and, then, its form has no con- 
straint. Usually, a simple ad-hoc term, ^V 2 u, is assumed 
and hence viscosity is considered as an artificial effect in 
the large-scale fluid, despite the good agreement of the 
model with numerical simulations (e.g. 0|). Actually, 
there are some alternative approaches (without introduc- 
ing viscosity) to the nonlinear regime like the linear po- 
tential and the frozen flow approximations (e.g. 
But a systematic comparison of several statistics for the 
distributions resulting from these approximations shows 
that the adhesion model works better over the nonlinear 
regime 6\. This situation suggests the need of a more 
physical motivation to invoke the viscous term. In this 
work, we use weakly nonlocal hydrodynamics to obtain a 
dynamical equation for the large-scale fluid in which the 
viscosity term is physically motivated. 



itj(r) = Ui(0) + r. 
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(1) 



Thus, 



Ui(r) - tii(0) « rj 



duj 



(2) 



The quantity = (Vu)y is called the velocity gradient 
tensor. The symmetric part of the velocity gradient is the 
rate of the strain tensor 



and the antisymmetric part is the vorticity tensor 

If dui duj\ 1 , . 

where the vector uj is the vorticity u> = V x u. Thus 



Uj(r) - Ui(0) k rjeij + -(w x r)j 



(5) 



We assume an irrotational large-scale astrophysical fluid, 
which is modelled as an incompressible Newtonian fluid 
with constitutive equation given by 
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where cr, 

shear tensor is linearly related to the velocity gradient 
tensor, we have, for an isotropic fluid 



< y is the stress tensor, p is the thermal pressure 

Assuming that the 
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with Aijki = XSijSki 
a symmetric tensor, 
comes 



+ ^ikhl + [J-'SuSjk and p = p! for 
Thus, the consitutive equation be- 



-pSij + 2fj,ei 



(8) 



in the right side of (13) is introduced just to prevent the 
inertial broadening of the first nonlinear structures - the 
Zeldovich pancakes. Hence, it is not a real viscosity , in 
the sense of a significant velocity gradient in the fluid, 
but a term which forces particles to stick together. 

Despite the use of a mock viscosity, the adhesion model 
shows a remarkable agreement with the large-scale struc- 
ture produced in N-body simulations 0- @ Thus, any 
simple improvement to this model does not need to 
change its mathematical structure very much, but just 
find a less artificial way to introduce the term in the 
right side of (13). Particularly, in this work, we want to 
find a physically motivated "viscosity" term presenting a 
dependence on the internal structure of the pancakes. 



Introducing this into the momentum equation 
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and making the spatial derivatives with respect to x = 
r/a(t), we find the Euler equation for the peculiar veloc- 
ity field 
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where pj p = v is the kinematic viscosity, <f> is the gravi- 
tational potential and, for the scales of interest, we have 
neglected the thermal pressure. 

Now, changing the time variable from t to b (the grow- 
ing mode of the linear theory) and rescaling the peculiar 
velocity and the potential gravitational as 
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ab 
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The correct framework to study continuum physics 
dealing with internal structures is the weakly nonlocal 
(coarse grained) hydrodynamics Q. In this context, we 
can use the Ginzburg-Landau equation as a first weakly 
nonlocal extension of a homogeneous relaxation equation 
of an internal variable. The Ginzburg-Landau equation 
is a modulational equation for a complex function, ip, 
which describes weakly nonlinear phenomena in contin- 
uous media with linear dispersion of general type. It is 
given by 
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ip+ (1 -Mc)V 2 V> - (1 + id)\4>\ 2 ip (14) 



where the parameters c and d are the linear and nonlinear 
dispersions. Some important properties of this equation 
deserve special attention. 

i) It can be reduced to a dissipative extension of the 
nonlinear Schrodinger equation 



we find 



zA7 2 v. 



(12) 



where fit) = 4 1" - , Qq = 8nGpo/3H§ and we have taken 
the velocity field as the gradient of a velocity potential 
(v = V#) as long as the motion is linear and streams of 
particles do not cross. 

Assuming that over the mildly nonlinear regime the 
gravitational potential is approximately equal to the ve- 
locity potential, we finally get 



i^ = VV±F(V#- 



(15) 



ii) It admits amplitude-phase representation in the 
form Q 



ip = i?exp i(. 



(16) 



iii) It can be derived from a variational principle. In 
particular, weakly nonlocal fluids respecting the second 
law of thermodynamics depends on the nonlocal potential 
given by 



— + v • V v = z/V 2 
ab 



(13) 



which is known as the adhesion approximation 0- For 
small v it reduces to the Zeldovich model which corre- 
sponds to a simple inertial motion. The viscosity term 



U = 5 ps{p,Vp) dV 



(17) 



where s is the entropy density |10| . 

iv) An important case, which allows us to connect hy- 
drodynamics and wave mechanics with a rich structure 
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of stationary solutions, is the Schrodinger-Madelung po- 
tential energy ffij 



U = 



ti 2 V 2 R 
2m R 



(18) 



IV. NONLOCAL ADHESION MODEL 



The wave mechanics analog of classical hydrodynam- 
ics is a well known subject in the literature (see e.g. [Tl| . 
[l2|). but the use of this equivalence as an approach to 
study structure formation in the universe is relatively 
new (e.g. [l^, 0)- The basic idea is very simple: 
to describe collisionless matter as a classical field ^{x, t) 
obeying the Schrodinger and Poisson equations 



- — V 2 + y)V, V 2 (/> = 47rGVV>*; (19) 
2m 



where V = m<f> and h is a parameter controlling the spa- 
tial resolution A, according to the de Broglie relation 
A = h/mv. This approach has the advantage of work- 
ing only with three spatial coordinates (plus time) and 
following multiple streams in the phase space. Promising 
results of this technique were found for the collapse of a 
self-gravitating object and to simulate a two-dimensional 
Cold Dark Matter universe [lflj . 

Here, we make the wave analog of a nonlocal fluid un- 
der a smooth gravitational field plus a coarse grained 
field by introducing the amplitude-phase ansatz tJj = 
Rexp(j-S) [15J into the Schrodinger equation (19). 
Equating real and imaginary parts separetely, we find 
two coupled equations for the classical action S and the 
real amplitude R: 



We can therefore attribute to the square of the amplitude 
of the wave function the interpretation of a fluid den- 
sity whose flux is conserved over time: p = tp*ip = \ip\ 2 - 
Then, we finally get 



(24) 



t()(x,t) = y/p(x,i)exp[iS/h]. 

with p — R 2 . 

Now, note that Equation (20) can be rewritten as 



du 

to— + m(u • V)u + V(V + U) = (25) 

or 



du 



(u- V)u + V( 



(26) 



where the quantity 



U 



v 

T 



V 2 p 1 fVp 
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(27) 



is the "viscosity" (or coarse grained) potential. Note that 
film has the dimensions of a kinematical viscosity, such 
that we have defined v = K/m. 

Using the same variable transformation we made in 
Section II for an expanding background, we finally come 
to 



^ + (v • V)v = A V (# - <p - ms ) (28) 

Again, for the weakly nonlinear regime, we have ip w 9, 
and then 
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(21) 



dv , _. 3 _ , 

^ + (W)v = --V^ 



(29) 



which is the weakly nonlocal adhesion approximation. It 
is important to note that the origin of viscosity should be 
related to the divergence of some "shear" tensor. Indeed, 
if we define 



where U is the potential previously defined in (18). Note 
that equation (21) can be rewritten as 



at to 



(22) 



kin 



T 



ViVj-p (v^xv.p) 



-VtVj lnp 



(30) 

we have that the divergence of this tensor is equal to the 
density force associated with 4>vis 



But p = VS in the Hamilton- Jacobi canonical transfor- 
mation. It then follows that u = VS/m and (22) takes 
the form of the hydrodynamical continuity equation 



V • a = — V0„i S . 
Now, defining a "kinematical velocity" as 



(31) 



dp 

dt 



V • (pu) = 0. 



(23) 



-Vlnp, 



(32) 
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the symmetric part of the velocity gradient (Vufcj„)jj is 
and 

thus showing the connection between coarse grained vis- 
cosity and nonlocal phenomena through the kinemati- 
cal shear tensor. The kinematical velocity was first in- 
troduced by 0] and it is not related to the organized 
mechanical motion, but to the perturbations due to the 
coarse grained nature of the fluid. From (34), we have 

— + (v-V)v = ^V 2 v fcm . (35) 

The kinematical velocity introduces a deviation with re- 
spect to the field v. If we define Av^n = v Vkin, we 
find the general expression 



— + (v • V)v = !/V 2 v + Vr? (36) 
do 

where Vrj = — i/V 2 (Avfcj„) is a noise term. Hence, the 
weakly nonlocal adhesion model corresponds to the noisy 
Burgers equation. Note that we recover the Zeldovich 
approximation if Vkin — > 0, which points out the impor- 
tance of internal structures to prevent inertial broaden- 
ing. 

V. DISCUSSION 



does not trace the internal dynamics of the pancakes. 
In this work, we assume the nonlocal hydrodynamics 
description in the amplitude-phase representation to in- 
troduce a physically motivated "viscosity" term into the 
adhesion model. Using the Schrodinger-Madelung poten- 
tial, we come to a new equation for the adhesion model as 
a function of a viscosity potential, which can be written 
as the divergence of a new shear tensor. 

Our model identifies the origin of viscosity as a nonlo- 
cal effect. Notice that the fluid equations are only valid 
when the average of the microscopic linear momentum 
equals the product of the averages of density and velocity. 
However, due to the mesoscopic nature of the averaging 
procedure, this is not guaranteed. Generally, local fluc- 
tuations in the fluid are small enough to fade away over 
averaging. But nonlocal (configurational) noise can add 
on the hydrodynamics noise, leading to non-selfaveraging 
properties in the matter flows. Mesoscopic arrangements 
in the system are sources of fluctuations, since the av- 
eraging procedure in the free volumes around them can 
lead to different results. This effect is a source of nonlocal 
noise ^(|. In the case of pancakes, the internal arrange- 
ments could be the spatial configuration of protogalactic 
gas clouds developed over the path followed by matter 
(inside pancakes) to form filaments and knots. These in- 
termediate structures can be considered as sources of the 
noisy mesoscopic fluctuations. 

Applications of (36) to the nonlinear regime and com- 
parisons to N-body simulations are beyond the scope of 
this paper and will be presented in a future contribution. 
For the moment, we conclude that the origin of viscos- 
ity in the adhesion model is probably associated with 
weakly nonlocal effects in the internal structure of the 
pancakes. These effects are given by some viscosity po- 
tential. In the case of significant mesoscopic fluctuations, 
we assume that the potential has a stochastic nature and 
the dynamics is given by the noisy Burgers equation. 



The wave-mechanical analog of hydrodynamics seems 
to be a promising framework to find alternative solutions 
to both analytical and numerical problems over the non- 
linear regime of structure formation. In this work, we 
have investigated the origin of the viscosity term in the 
adhesion model. Viscosity is the property which makes 
particles stick together once they enter into the caustics 
predicted by the Zeldovich approximation, just mimick- 
ing gravitational effects on smaller scales. This approx- 
imation describes the general nonlinear structures but 
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